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We present a search for the decay B + -> a^(1260)i : C* (892). The data, collected_with the BABAR 
detector at the SLAC National Accelerator Laboratory, represent 465 million BB pairs produced 
in e + e~ annihilation at the energy of the T(4S). We find no significant signal and set an upper 
limit at 90% confidence level on the product of branching fractions B(B + — ¥ a^(1260)A'* (892)) x 
S(a+(1260) ->■ 7r + 7r"7r+) of 1.8 x 10" 6 . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Measurements of the branching fractions and polar- 
izations of charmless hadronic B decays are useful tests 
of the standard model and a means to search for new 
physics effects. In decays of B mesons to a pair of spin- 
one mesons, the longitudinal polarization, fa, is partic- 
ularly interesting. Simple helicity arguments favor fa 
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to be close to 1, but several vector- vector (VV) decay 
modes such as B ->• <f>K* Q and B+ -> p+K*° % % 
are observed to favor fa ~ 0.5. Possible explanations for 
this discrepancy have been proposed within the standard 
model [4| as well as in new physics scenarios Q. 

New ways to explore the size of contributing ampli- 
tudes in charmless B meson decays and their helicity 
structure may come from measurements of the branch- 
ing fractions and polarization of charmless decays of B 
mesons to an axial-vector meson and a vector meson 
(AV) or to an axial- vector meson and a pseudo-scalar me- 
son (AP). Theoretical decay rates have been predicted 
with the naive factorization (NF) [6] and QCD factoriza- 
tion (QCDF) [7] approaches. The NF calculations find 
the decay rates of B — > AV modes to be smaller than 



the corresponding B — > AP modes. The more complex 
QCDF calculations find the reverse. For example, QCDF 
predicts a branching fraction of (HJ^ig q ) x 10~ 6 f° r 
B+ -> a+K*° and (32^;™°) x lO" 6 for B° -> &^ + , 
while NF predicts a branching fraction of 0.51 x 1CP 6 
and 1.6 x 10~ 6 , respectively. The first uncertainty on the 
QCDF prediction corresponds to the uncertainties due to 
the variation of Gegenbauer moments, decay constants, 
quark masses, form factors and a B meson wave func- 
tion parameter and the second uncertainty corresponds 
to the uncertainties due to the variation of penguin an- 
nihilation parameters. The NF prediction does not give 
an uncertainty on their value. 

B meson decays to charmless AV final states are sen- 
sitive to penguin annihilation contributions, which en- 
hance some decay modes while suppressing others. Thus, 
investigating decays to many final states will help deter- 
mine the size of the contributing amplitudes. 

A number of searches for AV decays to the final states 
at p~ , b\p and b\K* are presented in Ref. Q and Ref. 
|9j , with upper limits on the branching fractions of 30 x 
10~ 6 at 90% confidence level (C.L.) for afp~ and from 
1.4 to 8.0 x 10" 6 at 90% C.L. for the b x p and b x K* final 
states. In this paper we present a search for the decay 
B+ -> a+K*°. 

The data for this measurement were collected with the 
BABAR detector [1(| at the PEP-II asymmetric-energy 
e + e~ storage ring located at the SLAC National Acceler- 
ator Laboratory. An integrated luminosity of 424 fb~\ 
corresponding to (465 ±5) x 10 6 BB pairs, was produced 
in e + e _ annihilation at the T(45) resonance (center-of- 
mass energy y/s = 10.58 GeV). 

A detailed Monte Carlo program (MC) is used to sim- 
ulate the B meson production and decay sequences, and 
the detector response [ll|. Dedicated samples of MC 
events for the decay B + — > a^K*° with a^ — > p°n + 
and K*° — » K + ir~ were produced. For the af me- 
son parameters, we use the values given in Ref. [12| 
for studies with MC while for fits to the data we use 
a mass of 1229 MeV/c 2 and a width of 393 MeV/c 2 , which 
were extracted from B° — > a^ir^ decays [R||. We ac- 
count for the uncertainties of these resonance parameters 
in the determination of systematic uncertainties. The 
a^ — > 7r + 7r~7r + decay proceeds mainly through the inter- 
mediate states p°n + and an + [12|]. No attempt is made 
to separate contributions of the dominant P wave p° from 
the S wave a in the channel 7r + 7r~. The difference in ef- 
ficiency for the S wave and P wave cases is accounted for 
as a systematic uncertainty. 

We reconstruct af candidates through the decay se- 



The other primary 



quence af — > p u 7r + and p u 
daughter of the B meson is reconstructed as K*° — > 
K + tt~. Candidates for the charged kaons must have 
particle identification signatures consistent with those of 
kaons. Candidates for the charged pions must not be 
classified as protons, kaons, or electrons. We constrain 
the range of mass of reconstructed final-state candidates: 
between 0.55 and 1.0 GeV/c 2 for the p°, between 0.9 and 



1.8 GeV/c 2 for the a^, and between 0.8 and 1.0 GeV/c 2 
for the K*°. 

B + candidates are formed by combining af and K*° 
candidates. The five final decay tracks in a candidate 
are fit to a common vertex. Candidates which have 
a x 2 probability for the fit greater than 0.01 are re- 
tained. For these candidates, we calculate the energy 

substituted mass, tties = \/i s ~Pb7 an d the energy 

difference, AE — Eb — \\fs, where (E B ,p B ) is the 
B meson energy-momentum four-vector, all values be- 
ing expressed in the T(AS) rest frame. We keep can- 
didates with 5.25 GeV/c 2 < to es < 5.29 GeV/c 2 and 
\AE\ < 100 MeV. 

We also impose restrictions on the helicity-frame de- 
cay angle 9k*° of the K*° mesons. The helicity frame 
of a meson is defined as the rest frame of that meson, 
where the z axis is the direction along which the boost 
is performed from the parent's frame to this frame. For 
the decay K*° —> K + ir~ , 9x»a is the polar angle of the 
daughter kaon, and for a^ —¥ p°7r + , 9 + is the polar 

angle of the normal to the a^ — > 3ir decay plane. We 
define Hi = eos(#j), where i — (K*°,a^). Since many 
background candidates accumulate near |%^-*o| = 1, we 
require -0.98 < H K *o < 0.8. 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e + e~ — > qq events (q — 
u,d,s,c). We reduce this background source with a re- 
quirement on the angle 9t between the thrust axis [14| 
of the B + candidate in the Y(4S) frame and that of the 
charged tracks and neutral calorimeter clusters of the rest 
of the event. 

The distribution is sharply peaked near | cos #t | = 1 
for jet-like continuum events, and nearly uniform for B 
meson decays. Optimizing the ratio of the signal yield 
to its (background dominated) uncertainty, we require 
|cos6» T | < 0.8. 

A secondary source of background arises from b — > c 
transitions. We reduce this background by eliminating 
events in which one of the pions in the B + candidate is 
also part of a D candidate. 

Such D candidates, reconstructed from K~ir + and 
A'~7r + 7r + , are required to have an invariant mass within 
0.02 GeV/c 2 of the nominal D meson mass. 

The number of events which pass the selection is 15802. 
The average number of candidates found per event in 
the selected data sample is 1.5 (2.0 to 2.4 in signal MC 
depending on the polarization). 

We define a Neural Network (NN) for use in selecting 
the best B + candidate. The x 2 probability of the vertex 
fit and the p meson mass were the input variables to the 
NN. 

To further discriminate against qq background we con- 
struct a Fisher discriminant T [15j which is a function 
of four variables: the polar angles of the B + candidate 
momentum and of the B + thrust axis with respect to 
the beam axis in the T(AS) rest frame; and the zeroth 
(second) angular moment Lq (L2) of the energy flow, ex- 



eluding the B candidate, with respect to the B thrust 
axis. The moments are defined by Lj —J2iPi x \ C0S @i\ J > 
where Oi is the angle with respect to the B thrust axis of 
a track or neutral cluster i, and pi its momentum. 

We obtain yields and the longitudinal polarization /& 
from an extended maximum likelihood (ML) fit with 
the seven input observables AE, toes, ■?", the resonance 

masses m„+ and m K *n , and the helicity variables % K *a 

a i 

and H + ■ Since the correlation between the observables 
in the selected data and in MC signal events is small, 
we take the probability density function (PDF) for each 
event to be a product of the PDFs for the individual 
observables. Corrections for the effects of possible corre- 
lations, referred to as fit bias yield, are made on the basis 
of MC studies described below. The components in the 
ML fit used are: signal, qq background, charm BB back- 
ground, charmless BB background, and B + — » a^K* 
background. 

We determine the PDFs for the signal and BB back- 
ground components from fits to MC samples. We de- 
velop PDF parameterizations for the combinatorial back- 
ground with fits to the data from which the signal re- 
gion (5.26 GeV/c 2 < to es < 5.29 GeV/c 2 and \AE\ < 
60 MeV) has been excluded. 

For the signal, the toes an d AE distributions are 
parametrized as a sum of a Crystal-Ball function [16| 
and a Gaussian function. In the case of toes for qq 
and BB backgrounds we use the threshold function 
xy/l — x 2 exp [— £(1 — x 2 )\ , where the argument x = 
2toes/v / 'S and £ is a shape parameter. This function 
is discussed in more detail in Ref. 17]. In the case 
of AE for qq and BB backgrounds we use a polyno- 
mial function. The PDFs for the Fisher discriminant 
7-j(.F) are parametrized as a single Gaussian function 
or a sum of two such functions. The PDFs for the in- 
variant masses of the a+ and K*° mesons for all com- 
ponents are constructed as sums of a relativistic Brcit- 
Wigner function and a polynomial function. We use a 
joint PDF Vj{%k*°i^- a +) f° r the helicity distributions. 
The signal and the B + —> a^K* background compo- 
nent is parametrized as the product of the corresponding 
ideal angular distribution in Hk*o and H + from Ref. 118] 

times an empirical acceptance function <5(%jf*° >%„+)> 
while the helicity PDF for the other components is simply 
the product of the helicity PDFs for Hk*° and H a + . The 

Hi distributions for qq and BB backgrounds are based 
on Gaussian and polynomial functions. 
The likelihood function is 



N\ 



« 3 



V, (m* + )V 3 (to^.o )V 3 (H l K ,o ,Hl t ), 

where N is the number of events in the sample, and for 
each component j (signal, qq background, b — > c tran- 
sition BB background, charmless BB background, or 



TABLE I. Summary of results for B + -> a^ K*° . Signal yield 
Y, fit bias yield Yt., the branching fraction B = B(B + — > 
a^K* ) x 6(a + — > 7r + 7r - 7r + ), significance S (see text) and 
upper limit UL. The given uncertainties on fit yields are statis- 
tical only, while the uncertainties on the fit bias yield include 
the corresponding systematic uncertainties. The branching 



tion of K*° -> K+iT- 


is assumed to be |. 


Y Y b 


23(l(r b ) S UL (KT") 


61+2? 34 ±17 


0.7±gJ±?;| 0.5 1.8 



B + —t a^K* background), Yj is the yield of compo- 
nent j and Vj(x l ) is the probability for variable x of 
event i to belong to component j. We allow the most 
important parameters (first coefficient of the polynomial 
function for AE, the invariant masses of the af and the 
K* Q , and the width of the Breit-Wigner for the invariant 
mass of the K*°) for the determination of the combina- 
torial background PDFs to vary in the fit, along with the 
yields for the signal, qq background and b — > c transition 
BB background. 

We validate the fitting procedure by applying it to en- 
sembles of simulated experiments with the qq component 
drawn from the PDF, and with embedded known num- 
bers of signal and BB background events randomly ex- 
tracted from the fully simulated MC samples. By tuning 
the number of embedded events until the fit reproduces 
the yields found in the data, we find a positive bias yield 
Yt,, to be subtracted from the observed signal yield Y . 
The fit bias yield arises from the neglected correlations 
in signal and BB background events. 

The corresponding numbers are reported in Table [I] 

We do not find a significant signal thus we do not re- 
port a measurement on the quantity /&. In order to ob- 
tain the most conservative upper limit, we assume /l = 1 
in estimating the branching fraction. 

We compute the branching fraction by subtracting the 
fit bias yield from the measured yield and dividing the 
result by the number of produced BB pairs and by the 
product of the selection efficiency and the branching ra- 
tio for the B(K*° —} K + tt~) decay. We assume that 
the branching fractions of the T(45) to B+B~ and B°B° 
are equal, consistent with measurements [121 ] . The effi- 
ciency for longitudinally and transversely polarized sig- 
nal events, obtained from the MC signal model, is 12.9% 
and 18.6%, respectively. The results are given in Table 
HI along with the significance, S, computed as the square 
root of the difference between the value of —2 In £ (with 
additive systematic uncertainties included) for zero sig- 
nal and the value at its minimum. In Fig. [T]we show the 
projections of data with PDFs overlaid. The data plotted 
are subsamples enriched in signal with the requirement of 
a minimum value of the ratio of signal to total likelihood, 
computed without the plotted variable. We used 0.9 as 
requirement on the ratio in Fig. Q]for each variable. The 
efficiency of these requirements for signal is between 57% 




FIG. 1. Distributions for signal-enhanced subsets (see text) 

of the data projected onto the fit observables for the decay 

B + -> a+K*°- (a) m E s, (b) AE, (c) T, (d) m(p7r) for the a+ 

candidate, (e) m(KiT) for the K*° candidate, (f) T-L K »a and 

(g) W + . The solid lines represent the results of the fit, and 

a i 
the dot-dashed and dashed lines the signal and background 

contributions, respectively. These plots are made with re- 
quirements (see text) on the ratio of signal to total likelihood, 
computed without the plotted variable. 



and 70% depending on the variable. 

Systematic uncertainties on the branching fraction 
arise from the imperfect knowledge of the PDFs, BB 
backgrounds, fit bias yield, and efficiency. PDFs uncer- 
tainties not already accounted for by free parameters in 
the fit are estimated from varying the signal-PDF pa- 
rameters within their uncertainties. For K*° resonance 
parameters we use the uncertainties from Ref. [l2| | and for 
the a i resonance parameters from Ref. [13j |. The uncer- 
tainty from fit bias yield (Table QJ includes its statistical 
uncertainty from the simulated experiments, and half of 
the correction itself, added in quadrature. 

To determine the systematic uncertainty arising from 
our imperfect knowledge of the branching fractions of 
charmless B decays, we vary the charmless BB back- 
ground component yield by 100%. We conservatively as- 
sume that the branching ratio of B + — >• a^K*° could 
be as large as that of B + — > a^K*° and vary the 
B + —> a,2 K*° from to 18 events around a fixed yield 



of 9 events used for the B + — >• a^K* component in the 
likelihood function. 

The uncertainty associated with Jl is estimated by 
taking the difference in the measured branching fraction 
between the nominal fit (//, = 1) and the maximum and 
minimum values found in the scan along the range [0, 1]. 
We divide these values by y/3, motivated by our assump- 
tion of a flat prior for //, in its physical range. 

Uncertainties in our knowledge of the tracking effi- 
ciency are 0.24% per track in the B + candidate. The 
uncertainties in the efficiency from the event selection are 
below 0.6%. The systematic uncertainty on the measure- 
ment of the integrated luminosity is 1.1%. All systematic 
uncertainties on the branching fraction are summarized 
in Table HU 



TABLE II. Summary of systematic uncertainties of the deter- 
mination of the B + —¥ afK*° branching fraction. 



Source of systematic uncertainty 



Additive uncertainty (events) 
PDF parametrization 
a+ meson parametrization 
ML fit bias yield 
Non resonant charmless BB background 



B — > a%K charmless background 
Remaining charmless BB background 
Total additive (events) 



4 
6 

17 
3 
6 

7 
22 



Multiplicative uncertainty (%) 

Tracking efficiency 1.2 

Determination of the integrated luminosity 1.1 

MC statistics (signal efficiency) 0.6 
Differences in selection efficiency for af decay 3.3 

Particle identification (PID) 1.4 

Event shape restriction (cos St) 1.0 

Total multiplicative (%) 4.1 

+0.0 
-1.2 



Variation of f L [B(10 -6 )] 



Total systematic uncertainty [13(10 )] 



+0.6 
-1.3 



We obtain a central value for the product of branching 
fractions: 



B(B+ ->■ a+K*°) x B(a+ 



— > 7T + 7T 7T +N 



= (0.7 



+0.5+0.6\ 
-0.5-1.3/ 



x 10" 



where the first uncertainty quoted is statistical, the sec- 
ond systematic. Including systematic uncertainties, this 
result corresponds to an upper limit at 90% confidence 
level of 1.8 x 10~ 6 . 

Assuming S(aj t (1260) — > n + ir~Tr ± ) is equal to 
B{af (1260) -> n ± TT°ir ), and that B(af(1260) -> 3tt) 
is equal to 100%, we obtain a central value: 



B{B A 



aJK 



*o \ 



(1.3 



H.l+l.l> 
-1.0-2.6/ 



X 10" 



where the first uncertainty quoted is statistical, the sec- 
ond systematic. Including systematic uncertainties, this 



result corresponds to an upper limit at 90% confidence 
level of 3.6 x 10" 6 . 

This upper limit is in agreement with the prediction 
from naive factorization and lower than, but not incon- 
sistent with, that of QCD factorization. 
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